ATOMISTIC-TO-CONTINUUM COUPLING APPROXIMATION OF A 
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Abstract. We consider an atomistic model defined through an interaction field satisfying 
a variational principle, and can therefore be considered a toy model of (orbital free) density 
functional theory. We investigate atomistic-to-continuum coupling mechanisms for this atom- 
istic model, paying special attention to the dependence of the atomistic subproblem on the 
atomistic region boundary and the boundary conditions. We rigorously prove first-order error 
estimates for two related coupling mechanisms. 



1. Introduction 

The quasicontinuum (QC) method and, more generally, atomistic/continuum coupling (a/c) 
methods, are numerical coarse-graining techniques for the efficient simulation of phenomena 
and processes in materials at the nano-scale, such as defects, fracture, grain boundaries, or 
nano-indentation [T§\ [TE[ lllj . Incompatibilities between the treatment of forces in atom- 
istic and continuum models lead to difficulties in defining coupling mechanisms that do not 
introduce additional errors. Substantial effort has been made to understand this problem and 
to construct efficient and accurate a/c methods; see [T71 01 [151 El I2Q] for examples of formu- 
lations of computational methods and [TJ [2j [3J, [12| [131 [T^j and references therein for examples 
of analytical treatments. Formulations of a/c methods for atomistic models based on quantum 
mechanics were proposed in [U [6], but, to the best of our knowledge, no rigorous analysis of 
these methods exists. 

In the present article we formulate and analyze one-dimensional a/c methods for an atomistic 
model that is defined through an interaction field satisfying a linear variational principle. Our 
results are related to two classes of a/c methods: Firstly, our work can be viewed as an analysis 
of (a simplified version of) the a/c method proposed by Iyer and Gavini [9], who use field-based 
versions of classical potentials to formulate their method. Secondly, the atomistic model we 
formulate can be considered a toy model of (orbital free) density functional theory, and hence 
our work represents a preliminary step towards a rigorous analysis of the a/c methods described 
in [HE]. 

The article is structured as follows. In Section [T we formally motivate the atomistic model, 
and introduce the necessary notation. In Section [2 we give a precise formulation of the model 
with periodic boundary conditions and derive a "weak formulation" for the resulting forces on 
the particles. Section [3] is devoted to the analysis of the model in a bounded domain when 
the fields are subjected to Dirichlet boundary conditions. The Cauchy-Born continuum model 
is derived and analyzed in Section [4} Finally, in Sections [5] and [6] we propose two possibile 
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constructions of a/c methods based on different exchange of boundary conditions between an 
atomistic and continuum region, and establish error estimates. 



1.1. Field-based formulation of pair interactions. The following outline follows ideas 
presented in [9]. Let y = (yi, . . . , yjv) 6 ^ represent the coordinates of iV particles in one 
dimension. We consider an atomistic energy based on a pair-potential V, 

1 N 

The force on particle i is given by 

JV 

-D Vi S(y) = -^2sign( yi - yj )V' (\ yi - Vj \). 

i=i 

We note that the forces are nonlocal expressions in the sense that their computation involves 
summation over the other N — 1 particles. 

Next, we make a few modifications to this model. First, we replace the pointwise particles 
with smooth, nonnegative, and compactly supported particle densities 5 £ (- — yi) (such that 
J R 5 £ (x) dx = 1). This leads to 



1 ^ 

^(l/) ~ 2 Xi L l m 6 e(z - Vi)V(\z - x\)5 £ (x - yj)dzdx. 



To simplify the presentation further, we include the self-energies of the individual particle 
densities and define 



£ *(v) = o / / S £ (z-yi)V(\z-x\)5 £ (x-yj)dzdx. 

This additional self-energy contribution does not affect the forces. It can be computed explicitly 
and subtracted from the energy later on. Upon introducing the field (j) : K — > K, 

N 

<t>{x)= l p y (z)V(\x- z\)dz, where p y {z) = V) S e (z - yi), (1.1) 
to rewrite the energy £ £ (y) in the form 

£e{y) = « / Py(x)(f>(x)dx. 

z Jr 

It is now easy to see that the forces are given by the local expression 



-D y £ £ (y) = - / D y p y (z)(f)(z)dz. 
Jr 



Hence, if the field (ft is known, then it becomes unnecessary to compute nonlocal sums over 
particles. The nonlocality of the interaction has been encoded in the field (ft. However, it is 
now necessary to compute the field (ft, which is defined via the convolution (1.1). 

Suppose that the pair-potential V is the Green's function for a linear differential operator 
Ly(V); then, (ft can alternatively be computed by solving the differential equation 

L v (X7)(ft = p y . 
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As an example we consider the Yukawa potential in one space dimension 

FW = _L e -*i = _Ly K _2_ e ** dk . 

In this case can be obtained as the solution to 

— Acj) + m 2 (f) = p y 
or, equivalently, as a solution to the minimization problem 

(/> = argmin{- / \X7 (p\ 2 + m 2 f 2 dx — f p y ipdx 
- : I - ./ - Jr 

The resulting interaction potential £ e can also be written in the form 

£ £ (y) = -mm\l [ \V p\ 2 + m 2 <p 2 dx - [ p y ipdx\. (1.2) 
f 1 2 Jr Jr J 

The present work is devoted to the analysis of a/c approximations of (1.2) in a periodic one- 
dimensional setting. What distinguishes this analysis from previous analyses of a/c methods is 
that the coupling is achieved through an exchange of boundary conditions for the interaction 
field (f), rather than ghost-force removal ideas such as [H 115] . 



Remark 1. The interaction defined by (1.2) is purely repulsive. A purely attractive in- 
teraction can be obtained by changing the outer minus sign in the definition of £ e to a plus 
sign. We could combine two energies of the form (1.2) with different parameters m to model 
an interaction similar to the Morse potential = e _2 ' x ' — 2e _ ' x ' [9]. □ 

1.2. Notation. We consider an infinite chain of atoms on the one-dimensional lattice X = eZ, 
where e = 2/(2iV + 1) is the reference lattice spacing. Moreover, to keep the analysis simple, 
we admit only (2A^ + l)-periodic displacements from the reference lattice (cf. |14|). Hence, we 
define the spaces of admissible displacements and deformations, respectively, by 

U = {u G R z : u j+{2 N+i) = Uj Vj G Z, J2jL-N u j = °}> and 
y = FX+U, 

where F > is a prescribed macroscopic strain. A deformation y G y defines the computational 
domain 

^ = {V-n-uVn) 

for the field variable 4>. We note that the length of the interval is independent of y. 

We define the finite differences y' , y" G IA for y G y or U by their respective components 

/ _ Vj ~ // _ ~ 2 Vj + 

Vj £ > Vj £ 2 

Let us also define the weighted I 2 scalar product and norm by 

TV 

(•u 3 t>) e = £ ^2 u v v v Vu,v£U, \\u\\f2 := (u,u)l^ 2 Vit G U. (1.3) 

u=-N 

The ^°°-norm is defined in the obvious way 

IliilUoo = max \u u \ G U. 

u=-N,...,N 

The space IA equipped with the discrete Sobolev seminorm ||u||^i,a = will be denoted by 

IA 1 ' 2 and its topological dual space by IA~ 1,2 . The norm on U~ 1,2 is given by 

||T|| w -i,2 = sup Ti m • 
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y-N-i y-N y-N+i vn-2 vn-i vn 



Figure 2.1. Sketch of the basic atomistic problem: the field 4> is periodic in 
= (y_7v-i> 2/iv) and p y is a smooth particle density representing the atoms 
with positions given by y E y. 



For monotonically increasing y G y (which we will write as y' > 0) we denote by S(y) C 
H 1 ^) the space of continuous functions that are linear on every interval Qi = (yi~i,yi), 
i G {— N, . . . , N}. Furthermore, we define S # (y) = S(y) n H#(n) to be the subset of all 
periodic functions in S(y). 

2. Periodic Boundary Conditions 

We now put the field-based interaction potential that was outlined above in a precise math- 
ematical framework. Let the functional / : H^(Q) x y — > R be denned by 

I((p,y) = J ^e 2 \V(p\ 2 + ^m 2 ip 2 ^ dx — J p y (pdx, where 

p y (x) = ey^S^x - yj), and 5 e (x) = e~ l 5i{x/e). 

jez 

Here, <5i is a symmetric, nonnegative, regularized delta distribution with compact support 



[— where ?o > and ^8\dx = 1; see Figure 2.1 We will frequently refer to the 

paramter which is fixed throughout the paper. 

We then define the interaction potential £ : y — > M by 

£(y) = - min I(<p,y). (2.1) 



The respective minimizer (see Figure 2.1) 

(f> = arg min IUp, y) 

is the periodic solution to the Euler-Lagrange equation 

- e 2 A<p + m 2 4> = p y in 0. (2.2) 

Although <j> depends on y, we will usually suppress this in our notation. It will always be clear 
from the context, which configuration (j) belongs to. It follows from (2.2) and integration by 
parts that 

1 



J il 

To determine equilibrium configurations subject to a given external force / G U -1 ' 2 we need to 
minimize the total potential energy Ef : y — > IR defined by 

E f (y)=£(y) + (f,y) £ . (2.3) 

A minimizer y G y of (2.3) satisfies the following Euler-Lagrange equation in IA~ 1,2 : 

DE f (y) = D£(y) + f = 0. 



l 



In the following we analyze the derivatives of £. In particular, we obtain a "weak" formulation 
for the first derivative D£ that acts as a natural connection point for the coupling with a 
continuum model. 



Proposition 1.1. The potential £ : y — > R defined by (2.1) is twice continuously Frechet 



differentiable. The components of the first derivative are given by 



D yj £{y) = -e I V5 E (x - yj)4>(x)dx 
Jn 



(2.4) 



for j G {— N, ...,N — 1} and by 



D yN £{y) = -e / (W £ (x - y_ iV -i) + VS £ {x - y N ))(f>(x) dx. (2.5) 
Jn 



Proof. The proof of this result is standard and can be found in [S], for example. □ 

We stress the fact that the forces —D y £(y) are local expressions. To calculate the force on 
atom j it is necessary to know (j) in supp<5 £ (- — yj) but there is no need to sum over all remaining 
atoms. This nonlocality is encoded in the field 4>. 

Next we establish the weak formulation for the forces on particles. This very much resembles 
the structure of the continuum equations and will be the basis for the a/c coupling in Section 
[5j A version of this calculation was already shown in [7J, which used an interpolant for the 
displacement that is constant on the support of every 5 £ (- — yj). To avoid this restriction, we 
modify and extend the argument in [7J. 

For simplicity we assume that the supports of the densities of different particles do not 
intersect: 

supp<5 £ (- -yi) nsupp5 £ (- -yj) = Vi,j€Z, i / j. 

Since, |supp<5 £ (- — y«)| = e?o> this is equivalent to \yj — y%\ > £<?o for i ^ j or, if y is an increasing 
sequence, y'j > ?o for all j G Z. 



Lemma 1.2. Let y G y satisfy y' > % and let (j) 6 H^(Q) be the associated field, defined 
by (2.2). Let u = (uj)j^z S U be a test vector and u G S#(y) the periodic piecewise linear 
interpolant of u, that is, u(yj) = Uj for j G Z. Then, 

N 

D£(y)-u= D y .£(y)-Uj= / a y (x)Vu(x) dx, (2.6) 

j=-N ^ n 

where a y = o Vy \ + a y ^ and 

o- y ,i(x) = le 2 \Vcf>\ 2 - \m 2 (f) 2 + p y 4>, 

^ (2-7) 
a y>2 (x) = e 2^ <P(x)V5 £ {x - yj)(x - yj). 

j=-N-l 
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Proof. We begin by multiplying the derivative (2.4) for j £ {—N, . . . , N — 1} by the component 



D y .£(y)uj = - euj / V5 e (x - yj)<j>(x) dx 
Jn 



e I u(x)V5 e (x — yj)4>(x) dx + e / (u(x) — Uj)V5 £ {x — yj)4>(x) dx 
Jn Jn 



= e / 5 £ (x — yj)u(x)V(j)(x) dx + e / 5 £ (x — yj)<fi(x)Vu(x) dx 
Jn Jn 

+ e f (u(x) - Uj )V5 e (x - yj)(p(x) dx =: r[ j) + T 2 (j) + T 3 (j) . 
Jn 

Here we have used integration by parts but there are no boundary terms since u, 4> and p y 
are periodic on 0. Using (2.5) we obtain a similar expression for D yN £(y)uN- Summing over 
j = —N, . . . , N we obtain 

N 

DS{y)-u = D yj £(y)-u j =T 1 +T 2 + T 3 , (2.8) 

j=-N 

where Tj = Y2jL-N > * e {!> ^> From p y = £ J2jez ~~ ^ immediately follows that 

p y (x)i?!)(x)Vn(a;) dx. 

For Ti we can carry out the following rearrangements 

Ti= p y uV(j)dx= / (-e 2 A(/> + m 2 0)uV(/)dx 
'n Jn 

£ 2 V0A(/> + m 2 0V0)udx = - / V(-e 2 |V(/>| 2 + m 2 (j) 2 )udx 
n '2 Jn 



1 / (e 2 |V</>| 2 -m 2 2 )Vudx. 



Here, we have again used integration by parts and the periodicity of all functions involved. We 
deduce that 



2 - I 

n 



a y ^i{x)Vu{x) dx 



with o~y i as defined in (2.7). 

Before turning to T3 we first note that, since u is piecewise linear, 
X — U i 

u(x) = Uj H —( u j - u j-i) = u j + ( x — yj)^u(x) for x e Qj = (yj-i,Vj), 

Vj ~ Uj : 
X — U i 

u(x) = Uj H — - = uj + (x- yj)Vu(x) for x G Qj+i = (yj,y j+1 ) 

Vj+l Uj 



Hence, T3 in the above equation (2.8) can be written as 

N 

T 3 = £ ^ I <P( x )V S e(x ~ yj)(u(x) - Uj) dx 



j=-N-l " Q 

e ^2 j <t>{x)V8 e (x — yj)(x — yj)Vu(x) dx = e j a y ^{x)Vudx, 



j=-N-l 



n 



n 



with a y 2 as defined in (2.7), which concludes the proof. 



□ 



Remark 2. 1. In more than one space dimension the above calculations can be generalized if 
a triangular, respectively, tetrahedral mesh with the atomic positions as nodes is constructed. 
For example, this leads to 

°"y,iO) = {-h 2 \ V( t>\ 2 ~ 2 m V + Py<f) id + ^<t> ® V0. 



2. A closer look at the calculations in the proof of Lemma 1.2 shows that the weak form 
can be obtained for semilinear models — e 2 Acj) + F'(<p) = p y with any convex function F. 
Even a fourth-order model of the form e 4 A 2 (/> — e 2 Acf) + F'(<p) = p y admits a similar weak 
formulation. □ 

As already suggested in the introduction the Green's function for the differential operator 
— e 2 A + m 2 id acting on functions defined on R is given by 

= -Le-TM (2.9) 
2em 

We therefore get the following explicit formulas for the function values (j>(x) and V</>(x) for 

Proposition 2.1. Let y £ y and let (j> = arg min^^i ^ I(tp, y) be the corresponding 
interaction field. Then, for every x 6 $7, 

0(x) = / G £ (x - z)p y (z) dz = - 1 - T [ 5 £ (z - y k ) e'T^ dz, (2.10) 
Jr 2m t^ J ^ 

[ G e (x - z)Vp y (z) dz = - 1 - V f V5 £ (z - y k ) e"?!- 2 ! dz. (2.11) 



Proof. The proof of this proposition is similar to the one of [5j Thm. 2.1]; see also |10| Prop. 
2.4]. ' □ 

The following is a consequence of the simple exponential form of the Yukawa potential and 
some elementary properties of the exponential function in one dimension. Let yt, yj £ R satisfy 
Uj > Hi + £?o, so that the supports of particle densities representing the atoms i and j do not 
intersect. Then, 

/ / 5 £ (z - yj)e~T\ z ~ x \5 £ (x - yi)dxdz = I I 5 £ (z - yj)e~T^ z ~ x ^ 5 £ (x - y,i) dx dz 
Jr Jr Jr Jr 

= e -Tto-«*> / e-T^-y^5 e ( z - yj )dz- / e^te"*)^ - x ) dx 



= p 2 e-T(yj-y>\ (2.12) 

where we have defined 

f m r m r 

fi= 5 £ (x)e~£ x dx= / 5 £ (x)es x dx= / <5i(a;)e m3: dx. 
Jr Jr Jr 

Although we will frequently use this property, it is not essential for our reasoning. It merely 
makes some calculations more convenient. 
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a-L y-K 



FIGURE 3.1. The atomistic model in the domain O a with Dirichlet boundary 
conditions g = [g L g R ] T . 



3. Dirichlet Boundary Conditions 



In this section we consider a version of the model (2.1) in the domain O a = (o^or) C K 
subject to Dirichlet instead of periodic boundary conditions. This concept will be used later 
on, for the formulation of a/c methods, as the atomistic subproblem. We set a = [at ai?] T E R 2 
and Aa = cir — cll- Throughout Section [3] we think of y = (y-K, ■ ■ ■ ,Vk) as an ordered element 
of ft 2 1 such that < y-K < ' ' • < Vk < &R- The particle density p y is canonically defined 
by 

K 

Py = £ Yl 5 ^-Vj)- 

j=-K 

For simplicity we assume that the yj are separated and lie well inside O a in the sense that 
supply n dQ a = or, equivalently, 

y'i > %, for i = -K + 1,...,K, 
a-R-yK > e?o/2 and y- K - a L > e%/2. 



(3.1) 



We impose the following boundary conditions on the resulting field (j) : Q, a — > R: 

4>{o-l) = gi, 4>( a n) = gR; 

i.e., 4>\dn a = g with g = [g^ 5_r] t E M 2 . The interaction potential £ ag : Q 2K+1 — > M. is defined 
by 

£a, g (y) = ~ min I a (ip,y), (3.2) 
where J a : F 1 ^) x S1 2X+1 ^ R is given by 

r-a R ra R 

I a (f,y)= (\e 2 \V^\ 2 + \m 2 ^ 2 )dx- py^dx. (3.3) 

Ja L Ja L 

For given y the minimizer </> is the weak solution to 



— e 2 A(f) + m 2 (j) = p v in fi a , 



(3.4) 



<fe a = g- 

We will frequently use the decomposition 

4> = <f>0+£,a,g, (3.5) 

where (f>o E Fo(f2 a ) and £ aj9 E F 1 (r2 a ), respectively, solve the boundary- value problems 

-£ 2 A(p + m 2 (po = p y in Q a , 



and 



-e 2 ^ a ,g + m 2 ^ a>g = in Q a , 



(3.6) 



£,a,g\dQ a — 9- 

This last boundary-value problem can be solved explicitly, which yields the following lemma 



Lemma 2.1. The solution £ a;9 of (3.6) is given by 

Za, 9 ( x ) = c L (a,g)e-^ ( - x - aL) + c R (a,g)e- 
where the coefficients CL(a,g) and c^(a, g) are given by 

c{a,g) 

and we have defined r = exp(— ^Aa) 



c L (a 


9) 




"l T 


-1 


9L 


Ci?(a 


9). 




T 1 




9R 



l -(a R -x) 



-■■T^-g 



(3.7) 
(3.8) 



Note that, for Ao 3> e, r is exponentially small; hence we will often neglect terms of that order 
of magnitude. We will write 0(t) for a quantity or function that is (uniformly) bounded above 
by Ct in modulus, where C is independent of e and Aa. For example, we have c{a, g) = g+0{r). 

Next, we compute the derivative of £ a ^ g with respect to the atomic coordinates. For these 
derivatives, we obtain a "weak" formulation of the same shape as in the periodic case (see 
Proposition 

If y' > 0, then we denote by S(y L) a) the set of continuous, piecewise affine functions over 
the mesh given by the nodes a^, y-K, ■ ■ ■ , UK, clr- Moreover, So(y U a) = S(y U a) n Hq (fi a ). 



Proposition 2.2. Let a,g£ M 2 , < an; then E a ^ g : y 
Frechet differentiable at y. 

(i) The components of the first derivative are given by 



defined by (3.2) is continuously 



Dyj£a, g (y) = —e / V5 £ (x - yj)4>{x) dx for i = -K, ... ,K. 



(3.9) 



(ii) Let u €U be a test vector, u £ So(y U a) its interpolant, and let rainy' > then 



D y £ a ,g(y)-u 



a y {x)Vu{x) dx, 



(3.10) 



where o~ y is given by (2.7). 



Proof. The derivatives with respect to the coordinates y are easy to calculate along the same 
lines as in the proof of Proposition [Ll} The weak formulation can be obtained as in the periodic 
case (Lemma 1.2) using the fact that the interpolant u vanishes on d£l a . □ 



Remark 3. We point out that, in general, 

£a,g(y) + \ I Py<t>dx. 
z Jn a 

However, we will see below that £ a ^ g {y) can be written as the sum of a boundary data contri- 
bution and a term that is independent of g. □ 

With a view to the subsequent derivation of a/c methods we will from now on interpret a 
and g as arguments to £ a ,g rather than fixed parameters entering its definition. We consider 
the map x M. 2 x M. 2 — > M., (y, a, g) i— > £ a ^ g (y), and derive the derivatives of this map with 

respect to the boundary a and the boundary data g. 
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3.1. Dependence on the boundary positions. When formulating a/c methods in Section 
[5] we will let the boundary a of the atomistic subdomain depend on the configuration y. It 
is therefore necessary to understand the dependence of the energy £ a , g (y) on a. Our main 
result is that the derivative D a £ a ^(y) can be combined with D y £ a>g (y) into a weak formulation 



reminiscent of (2.6). This will be a central building block for a/c methods. 



Proposition 3.1. Suppose that y G y , miny' > Let h = [hi h R ] T G M. 2 and u = 
(u_X) • • • , uk) G M 2A+1 be test vectors, and let u G S(y U a) denote the interpolant of u and h 
in the sense that 

u(a L ) = h L , u{a R ) = h R , and u(y,j) = Uj Vj G {-K, . . . , K }. 

Then, 

D a £ a ^ g (y)-h + Dy£ a , g (y)-u = \ a y (x)Vu(x) dx. 



Proof. This is a direct consequence of the following two lemmas. □ 

In the first auxiliary lemma we compute the derivative of £ a ,g{y) with respect to a = [ol o,r] T 
while keeping the relative distances between the atoms constant. In other words we consider the 
change in £ a ^ g {y) when the whole domain Q, a is stretched with the atom positions following this 
stretching. For y G £l 2 a K+1 let X = {X- K , X K ) G (0, l) 2K+1 be given by yj = a L + AaXj 
for all j G {—K, ■ ■ ■ , K}- For fixed g, X we define 

£(a) := £ a , g (aL + («R - o>£)X), and 

D aR £a,g{y) ■= D aR £{a). 

(We understand cll + (or — ai)X in a componentwise manner: (<zl + AaX)j = ol + AaXj for 
all j G {—K, ... , K}.) The derivative D aL £ a ^ g {y) is defined analogously. 

Lemma 3.2. Let y G tt 2 a K+1 satisfy ( [371] ); then 

- D aL £ a)g (y) = D aR £ ai g(y) = -r- / 0-3,(2;) dx. (3.11) 



Proof. We fix X and let 77(a) := at + AaX. We begin by transforming the problem to the 
unit interval (0, 1) using the transformation x h-> X(x) = (x — a£)/{p,R — ol): 



£{a) =£ a , g (r 1 {a)) 



\e 2 W4 



\m 2 4> 2 + p v{a) 



dx 



Aa 



2Aa- 



-|V0p 



m 



+ Pr)(a) 



dX. 



1.1 



(3.12) 



that, to 



Here, 4>{X) = 4>(x(X)) and p^u^X) = p v ^(x(X)). It follows as in Proposition 
compute D a £(a), it is sufficient to calculate the partial derivatives of the right-hand side with 
respect to a R (the derivative of 4> or <f> with respect to a R does not appear since </> is a minimizer 
of L a (-,y)). This leads to 



D aR £(a) 



£ 



2Aa 2 
1 



2 

m 22 



+ Pr,{a)<l> )dX + Aa 



1 e 2 



Aa 3 



\V0\ 2 dX 



+ Aa I ^D aR p^ a) dX. 
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Transforming the first two integrals on the right-hand side back to the interval Q a we arrive at 

.2 

-P.. JiA 4- - 

Ac 



1 e 2 f 1 f 

\a Aa Jo Aa J n 



where a V) i was given in (2.7). 

It remains to differentiate p v (a) with respect to a R . By the definition of the transformation 
x i — y X(x) we have 



K 



Da R pr,(a)(X) = eD aR ^ 5 £ ((ttjj ;- Cl L ) (X - Xj )) 



K 



j=-K 



Using Aa(X — Xj) = (x — yj) we therefore get 



Aa $Da B p v ra)dX = -r- Y] / (x - yj)V5 e (x - yj)<f>(x) dx 

JO Aa ~~K Jn " 



a / ^y,2(x)dx 



with <Jy.2{x) as given in (2.7). 

To see that D aL £ = —D aR £ we simply note that £ (a) depends only on Aa, which can be 
seen from (3.12) and the definition of "p v f a \(X). □ 



We define 8 R £ S(y L) a) to be the piecewise linear function with 

9 R (a R ) = l, 8 R (a L ) = 0, e R ( yj ) = for alii G {-K,...,K}. 
The function 8l G S(y U a) is defined analogously. 



Lemma 3.3. Let y G satisfy (3.1); i/ien, i/ie derivatives of £ a , g (y) with respect to a^, 

a R (f or fixed y and g ) satisfy 

D aL £a,g{y)= I CT y (x)V8 L (x)dx, 
D aR £a,g{y)= \ <?y (?) V8 R (x) dx. 



Proof. Let Q R be the affin e function defined on fl a with Sij(ax) = 0, @ R (a R ) 

yields 



VQ R (x) = Lemma 



3.2 



D aR £a )V (y)= I o-yV@ R dx 



<jyV(G R -8 R )dx+ / a y V8 R dx. 

r2 a J S7 a 



Now, we have Q R — 8 R G So (y U a) and hence, by Proposition 2.2 

A" 



/ c ry (x)V(6 R -^)drc= D Vj £ a!g (y)Q R (y 

JQ-a A— IS 



1. Since 



(3.13) 



(3.14) 



j=-A- 
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However, D aR £ a g (y) was defined as the derivative with respect to a R , while the relative dis- 
tances of the atoms are kept constant. This can be formulated as 



A' 



D aR 8 a) g{y) = D aR £ a , g {y) + ^ D y .£ a ,g(y)@ R (yj). 

j=-K 



Inserting this into (3.13) and using ( |3. 14 ) then gives 

/ a y (x)V9 R dx = D aR £ ajg (y). 
Similarly, we can show the expression stated for D aL £ a ^ g (y). 



□ 



3.2. Dependence on the boundary conditions. Next, we compute the derivative of £ a ^ g (y) 
with respect to the boundary conditions g when the configuration y and the boundary a are 
kept fixed. We define 

j L (y,a) = 2 p y (x)G £ (x - a L ) dx, and j R (y,a) = 2 p y (x)G £ (a R - x) dx. (3.15) 

J S2 a J Vl a 



Lemma 3.4. The partial derivative of £ ay g{y) with respect to g is given by: 

. T 



Dg£ a ,g(y) = -mel (1 - r 



Ch{a,g) 
c R (a,g) 



lL(y,a) 



■ T~ 



where T a , c(a,g) = [ci{a,g) c R {a,g)Y and r = e e Aa are defined in Lemma 



2.1 



Proof. Throughout the proof we suppress the arguments of 7x, 7^, and c for ease of readability. 
We recall the additive decomposition <f> = 0o + £, a ,g from ( |3.5[ ). From <pQ G Hq(Q), and from 
the equation — e 2 A£ aj9 + m 2 ^ a ,g = it follows that £ 2 (V£a i9 , V<po) + fn 2 (£ Q)9 , <fto) = 0. Hence, a 
short calculation shows that the energy £ at g(y) can be rewritten as 



£ a ,g(y) = - I a{4>,y) = -i a {4>o,y) - h{U,g,y) 



(3.16) 



The first term on the right-hand side does not depend on the boundary conditions g and the 
second term is known explicitly: using — e 2 A^ a>9 + m 2 £ ajg = 0, integration by parts, and the 
explicit formula (3.7) for £ a „, we obtain 



I a (€a,g,y)= l(e 2 \S7£a, g \ 2 + rn 2 (£ tg )dx- / p y ia,gdx 



-£,a,g{a L )X7£ a ,g(aL) + £a,g{a R )V £a,g(a R )) 



em 



- me 



- me 



2 1 2\/i -2 — Aa\ 
C L + C R) I 1 - e E ) 



a, 



Py{c L e 



l -(x-a L ) 



Py^a,g dx 

+ c R e-T( a R-*))dx 



{l-r 2 ) 



C L + C R (1 _2\ 
2 

4 + 4 



2me 



Py [c L e 



+ c R e-f^-^)dx^ 



[l-T 2 ) - (c L J L + C R J R )^j . 



Here we have used the Green's function G £ from ( |2.9| ). Differentiating this expression with 
respect to cl and c R and applying the chain rule with D g c = T" 1 yield the result. □ 
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Remark 4- 1. We remark that D g £ a ^ g (y) = if and only if cr,(a, g) = 7i(y, o)/(l — r 2 ) and 
cn(a,g) = 7fl(y,a)/(l — r 2 ). According to (3.8) this corresponds to the "optimal" boundary 



conditions 



1 7 L + T7ij „ 1 T7L + Jr 

9l= i ~ i T ~ » and £r = 1 ~ -I I (3-17) 

1 — Tl + T 1 — Tl + T 



That is, the boundary conditions are weighted averages of the values jh^lL and j-~^jr. 



2. As can be seen from Lemma 3.4 the boundary data contribution I a (£a,g,y) to the energy 
£a,g(y) is quadratic in g. For fixed configuration y and domain £l a the boundary conditions 
g = g*(y,a) minimize the boundary data contribution I a {£a,g,y) to the energy £ a , g (y)- This 
is equivalent to minimizing I a (-,y) over iJ 1 (fi a ) and therefore leads to homogeneous Neumann 
boundary conditions for on dVt a . 

3. If Aa 3> e, i.e., t<1, then we have 7i/ij = + 0( T )> and hence we can simplify 

I a (Za,g,y) = me{\{gl + g%) - {g L 9L + 9R9* R )) +0(er), and (3.18) 
D g £ a ,g(y) = me(g* - g) + O(er). □ 

A useful auxiliary result for the analysis of a/c methods is the global Lipschitz continuity of 
the field <ft with respect to variations in the boundary conditions g. 

Lemma 4.1. Let 4>\,<j)2 G H l {Q, a ) be minimizers of I a (-,y) subject to the boundary conditions 
gi G M 2 , respectively, g2 £t 2 . Then, 

|0i(x)-0 2 (x)| < v^lT-^-^le-e^, and 

e|V0i(x) - V0 2 (x)| < IT- 1 ^ - ff2 )| e~ 

where d a (x) := min(x — ai,aR — a;) denotes the distance to the boundary of£l a , for x G fi a . 

Proof. We write both functions in the form ^ = </>o+£a,gi, « G {1)2}. For z = 1, 2, let q = T~ x gi 
be the respective coefficients entering £ a . 9i ; then 

|0iO*O- 02O)| = |^a, 9l (x)-ea, 92 0)| < |c ll i-c 2l L|e-T(™) + | Clifl _p2,H|e-T( «- !B ). 
This immediately yields the first bound. The bound for the derivatives is obtained similarly. □ 



3.3. A special case. We now take a closer look at the interaction potential £ a ^ g from (3.2) 
with the y-dependent boundary conditions g = g*(y,a) defined in Remark [4j 

Proposition 4.2. Let y G ft 2X+1 . Then, 

£a,g*(y,a)(y) = I J Py 0)e~ T Py (z) dz dx + tM t (j L , Jr) 

+ / Py( x )( e ~^ (2aR ~ X ~ Z) +e-f (x+z - 2aL) ) Py (z)dzdx, 



where M t {^i^r) depends quadratically on 7^ and 7^. 



Expression ( 3.19| ) can be interpreted as the energy of the atoms represented by y interacting 



with each other plus the interaction with mirror atoms outside Q. a . This mirror interaction was 
introduced by means of the boundary conditions g = g* . 

For the proof of the proposition it is convenient to use an explicit formula for the function 



values of <po G Hg(f2 ) from the decomposition (3.5). By Proposition 2.1, the Green's function 



for the equation — e 2 A0 + m 2 (p = p y in E is given by G £ (x,y) = e ' a y '- We will now 
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construct the Green's function G £>a for the operator — e 2 A + m 2 id subject to homogeneous 
Dirichlet conditions on dQ, a . 



Lemma 4.3. Let 0o £ ^o(^a) satisfy — e 2 A0o + m 2 (j)Q = p y in Q a . Then, 

00 (jc) = / G £5a (x,z) / Oy(2;)d2; Vx G O a , 

w/iere G £i a = Gri + rG^, urai/j G £ %, i = 1, 2, gwen fry 



2rae 



1 1 

2me 1 — t 2 



Te -- £ -( x + z ~ 2a L) _|_ Te -~r( 2a - R ~ :c ~ 2 "' 



7/t , V III f , * 



(3.20) 



Proof. The proof of this result is standard [5J Chapter 2.2.4]; see also [10} Lemma 3.10]. □ 
We remark that G £ , a = gQ + 0(t). 



Proof of Proposition ^ -£\ We have already seen in (3.16) that for any choice of boundary data 
g G M 2 the energy £ a ^ g (y) can be written as the sum of two terms 

£ a ,g(y) = -i a {4>,y) = -4(0o, y) -l a (€a, g ,y), 

where I a (<fio,y) is independent of the boundary conditions. 

Calculation of I a {(f)o,y) ■ Since the function 0o is a minimizer of I a (-, y) over Hq(Q), we have 
with the expression ( 3.20| ) for <fto(x) that 

Iff . . 1 



i a (<t>o,y) 



p y 0o dx 



f2 Q 



p y (x)G £) a{x, z)py(z) dz dx. 



(3.21) 



By the definition (3.15) of 7^ and 7^ we have 



1 



4me 



4me 



Pi/ We e 



(x+z—2a L ) 



p y (z) dx dz 



me 



"7l, 



(3.22) 



p y (x)e-T( z -^-^) /93/ (z)dxd Z = — 7L7R- 



Inserting the expression G £ , a = Gi]l + tG £ % into 03.211) and using these equalities yields 

1 



,(2) 



4(00,2/) 



p y (x)G E (x, z)py(z) dz dx + ^ ( 7 f, + 7!) 
(t7l + T7fl - 27L7fi) • 



+ 



me r 



Calculation of Ia(£ a ,g*(y,a)>y)- Prom Lemma 3.4 we know that for general g G 



/c 2 + c 2 \ 
I a {ta,g,v)=mel — - fl (l - t 2 ) - (c l7l + c r7r ) J . 

If (7 = g*(y,a), then cl = 7l/(1 — t 2 ) and cr = 7_r/(1 — t 2 ) as seen in Remark[4j Hence, 

me 1 



4(£a,g*(y,a), 2/) 



2 l-r : 
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Py 



Vj-i 



Qj+2 



26+2 



Figure 4.1. The Cauchy-Born approximation: independent periodic problems 
are solved on the cells Qj = (yj-x, Vj) leading to locally defined fields ifrV). 



Isolating the dependence on r gives 

UL^y^y) = -™tf L + - ^^(7l + 7l). (3-23) 
Conclusion. Adding —I a {(,a,g*(y,a),y) as J us t obtained and —I a ((j)o,y) from above we arrive 

at 

Iff -™\ x - z \ me , o o\ 

£a,g*(y,a)(y) = J— / Py(x)e ^ Py (z) dz dx + — ( 7i + J R ) 

(3 24) 

(T7L + IlLlR + T7r) . 



S7 a J 

me t 



4 1-t 2 



Defining tM t (jl,jr) to be the third term on the right-hand side and applying (3.22) yields 
pTT9| ). □ 

4. The Cauchy-Born Approximation 



The next building block for the design of a/c methods based on the model (2.1) is the 
respective continuum model. Let y € y satisfy mmy' > The Cauchy-Born approximation is 
obtained by computing the energy of the cells Qj = (yj-i,yj) independently from one another, 
by treating each of them as part of a homogeneous chain (see Figure 4.1). We define the 
Cauchy-Born energy of the cell Qj by 

Sf(y) = - min I I Ue 2 \Vip\ 2 + ^m 2 ^ 2 ) dx - [ p v ipdx). (4.1) 
^hI{Qj)\Jq 3 z Jq, J 

Note that this energy only depends on the distance (yj — yj-i)- The minimizer ijj^ of (4.1) 
satisfies the equation — e 2 A^^ + m 2 ip^ = p y in Qj and its \Qj\ -periodic extension to R: 

- e 2 Aip^ + m 2 ip {j) = p y(j) in R. (4.2) 
Here we have defined the positions yd) = (y^)kez of an infinite chain of equidistant atoms by 

y? = %' + (^-i)(%'-%'-i) vfcez. (4.3) 

The Cauchy-Born approximation £ ch (y) of the atomistic energy £(y) is then given by the sum 
over all cells 

N N 

£ ch (y) = £?(y) = 2 E / Pv*®**. (4-4) 

j=-N j=~N J Qi 

In the Cauchy-Born model we seek to minimize the total potential energy E c ^ : y — > R defined 
by 

Ef(y) = £ cb (y) + (f,y) £ . (4.5) 
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Whether the Cauchy-Born model is a good approximation to the exact atomistic model strongly 



depends on the regularity properties of minimizers of (4.5) 



Let u G U be a test vector and u G S#{y) an interpolant of u, i.e., u{yj) = Uj for j G Z. It 
follows as in Lemma 



3.2 



that the derivative of £f(y) can be written in the form 



D v £f(y)-u= U [~_^ 1 I of y (x)dx= I af y (x)Vu(x) dx, 



Uj — Uj-l 



Vj ~ Vj-l JQj 

where the local continuum stress function o~f y , in direct correspondence with (2.7), is 



(4.6) 



of Or) = ±e z \V^)(x)\ 2 - \m l ^\xf + p y (x)^>(x) 



U), 



N 



+ e E ^ (j) {x)V5 £ {x- yj ){x- yj ). 

j=-N-l 



(4.7) 



Furthermore, we define the Cauchy-Born stress function o y : Cl — > R by 

of{ x ) = of y {x) if x G Qj 

for all x G CI. 

4.1. Consistency. Next, we turn to the consistency analysis of the Cauchy-Born approxima- 



tion, for which we thoroughly analyze the modelling error incurred. From (2.6) and (4.6) we 
deduce that 

\(D£(y)-D£ cb (y))-u\< [ \o y (x) - o y b (x)\ \Vu(x)\ dx 

Jn 

N 

E J M X ) - a fy( x )\Nu(x)\dx, 



where the stress functions a y and af y are given by (2.7) and (4.7), respectively. To investigate 
the modelling error |<7 y (2;) — £r?y(a:)| incurred by going from the atomistic description to the 

Cauchy-Born approximation it is therefore sufficient to analyze \4> — t/j^'\ and |Vc/> — V?/>^| in 
Qj for every j G {— N, . . . , N}. 

Lemma 4.1. Let y G £°°(Z) and define yC?) = (y^)kez by y^ 1 = yj + ey'j(k — j) for all k G Z; 
then 

\y n ~ \ <(n- j)e?\\y"\\p(\j,n-i]) for n > j, 

\y n ~ Vn ] \ < (j - 1 - n)e 2 \\y"\\ i i {[n+l j_iY ) , for n < j - 1. 

Proof. Assume, without loss of generality that n > j. Since yj—\ = V^li an d yj = Vj*\ 

n n k— 1 n k—1 

yn- y $=e e {v'k-^ ) )') = e 2 E EU'-^rw E E^ 

k=j+l k=j+l l=j k=j+l l=j 

where have used that (y^')' is constant. Changing the order of summation we get 

n— 1 n n—1 

\y n -y^\<e 2 Y E \y'i\ = < (n-j)£ 2 \\y"\\m,n-i)y □ 



l=j k=l+l 



1=3 
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In the next result we estimate the errors \<j>(x) — ij}^(x)\, |V^(x) — V^^'^(a;)| for x in the 
cell Qj. As anticipated by Lemma 4.1 they depend on the second difference y" . 

Lemma 4.2. Let y G y satisfy miny' > Let (f> G fli(fl) satisfy (2.2) and tp^> G H\(Qj) 



satisfy (4.2), respectively. Then, 
|U- 



< 



ne- mnminy , and 



n=l 
oo 



n=l 



Proof. From Proposition 2.1 we immediately deduce that, for all x £ Qj, 



(4.8) 



dz. 



(i) (i) 

Since = yj and = y^-i, the respective terms in the sums cancel. Hence, we get for 
x G Qj: 

4>{x) - ^\x) = J2 f {5 £ (z - y k ) - 5 £ (z - !,£>))e-Tl*-*l dz. 

Zm /TED 



feez 



We now derive bounds on the individual terms in the sum. Note that (2.12) simplifies the 
following calculations but due to the smoothness of the Green's function similar bounds can be 
obtained without it. 



Let k > j. Then we have \x — z\ = z — x for all z G supp<5 £ (- — yk) and all z G supp5 e (- — y k ~ 
Thus, with ( |2.12[ ), 



1 

2m 



{8 e (z-y k )-8 e (z-y®)) 



e e 



|x-z| 



dz 



2m 



e-lW-e-TW -*>). (4.9) 



If J/* > J/fc, then 

1 

2m 



(<5 e (z- yfe )-«5 £ (z-^))) 



e e 



dz 



< e e 

2m 



At -m, y x) m (i) 
< — e e VM ' — (y) —ilk)- 
~ 2m e K!Jk y ' 



Using {yk — x) > (k — j)emmy' for all x G Qj and applying Lemma 4.1 leads to 

2/Jfe - V k 



_fl e -y(3/fe-a;)| 

2e 



■ Wl <^||l/ , lh tt »-i])(fc-j)e- (fc - i)mn,tay '. 



The same bound on (4.9) can be obtained if y k < y k - 

For any k < j — 1 we can use the same techniques to obtain that 

^ / {8 £ {z-y k )-8 £ {z-yf)) e -?l~l dz 



< 



12/ 



i J -i])(j-^-l)e-°- fe - 1)mmin3/ '. 
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Summing over all k G — 1, j} we deduce that 

oo 

\^{x)-^\x)\ < »eJ2\\y%H\j-n, j+ n-i])ne- mn ™y'. 

n=l 

The proof for the derivatives V</>, Vtp^ is analogous. □ 

We wish to prove modelling error estimates on — o*j y||nx>(Q .) i n terms of ||</>— || Loo( -q ^ 
and ||V<^> — ^^(q.)- Since the stress functions a v and a cb y are quadratic in the fields <p 

and ifju) we need L°°-bounds, which we establish in the next lemma. 

Lemma 4.3. Let y 6 y, y' > a^rf Zei <^> = argmin^g^i /q\ I((p, y) be the corresponding 

field. Then, there are continuous functions Kq, K\, that depend implicitly on m (but are 
independent of e and y ), such that 

||0lk°°(n) < Ko(mminy'), and e\\V<t>\\ L ^(a) < A'i(mminy')- 



Proof. The stated estimates follow in a straightforward manner from the integral representation 
of the solution <f>; see \10\ Lemma 4.4] for the details. □ 

We can now prove the following modelling error estimates. 



Lemma 4.4. Let o~ y and a cb y be given by (2.7), respectively, (4.7); then 



where the constant C only depends on 5\, K l = Ki{my'), and on m. 

Proof. From the definitions of the atomistic and continuum stress function we deduce that 

a y (x) - <r£ y (x) = - l{eV(/)(x) - eVV> (i) (x)) (eV0(x) + eV^ (j) (x)) 

+ \m 2 {ct>{x) - ip (j \x)) (<f>(x) + ip ij \x)) - p y {x)((j)(x) - 7p( j \x)) 



(4>(x) - ^{x)) ^ eV ^( S - Vi)( X ~ Vi 
i=j-l 



for all x £ Qj. With 5 £ (x) = e 1 5i(x/e), the L°°-bound on <f> from Lemma 4.3, and the 
analogous bound for we get 

!|eV0(z) +eV^^(ar)| < Ki(mmmy'), 
^-\4>(x) +tp^\x)\ < m 2 Ko(mmmy'), 
\\Pv\\l°° < ll^illi 00 ) 
\eV8 E (x-yi)(x-yi)\ < || VJiid||z,oo, 

which implies the stated result. □ 



is 



4.2. Stability. Besides consistency, the second crucial property of an approximation to a given 
model is its stability. The following auxiliary result will play a role in the stability analysis of 
a/c methods. 

Lemma 4.5. Let y G y satisfy miny' > ?o- Then, for all j G {— N, . . . , N}, 

2 2 

D 2 Sf (y)-[u,u\ > r ^J^e- m ™ xy 'e\u , j \ 2 Vw G W. 



Proof. We first recall that £j b (y) = \ Jq. pyip^ dx because ip^' is a minimizer of (4.1). Ex- 
tending ifrV' | Qj | -periodically to M and using the symmetry of the cell problem, we can rewrite 
this as 

2 Jr 

We now insert the explicit formula ( |4.8[ ) for tp^\x) and apply (2.12) to get 
£?(y) = ^T,f I Se(.x- yj )S e (z-y^)e-T\'-'\dz6x 



4m 2m 

fcez u=l 

where the constant £ se if coming from k = j in the sum represents the self-energies of the atoms 
in the cell Qj. Here we have also used that \y^ — yj\ = \k — j\y'j for all k G Z. Differentiating 
twice leads to 

2 00 

D 2 £f ).[u,u] = rn ^eY,^~ Vmv ' i W j ? 



u=l 

2 00 2 

^ m//^ , ; | 2 ^2 e -i/mmaxt/' ^ g — mmaxy' i / |2 

— 2 J ^ — 2 J 

!/=l 

In the last step we have only kept the term for v = 1, which represents the nearest neighbour 
interactions. □ 

5. Atomistic-to-Continuum Coupling 
The computation of the original atomistic energy £{y) involves the solution of the optimiza- 



tion problem (2.1) posed in the whole of = (y-N-i, Vn)- Our goal is the construction of 
computationally cheaper, approximate energies £ ac (y) such that £(y) ~ £ ac (y) for all relevant 
y and minimizers ?/ ac G y of 

Ef{y)=£^{y) + (/,y) ej 

are good approximations of minimizers y of the energy Ef from ( |2.3[ ). 

Following the philosophy of a/c methods we approximate £(y) by the continuum model 
where y is smooth and a version of the atomistic model where y is nonsmooth. In the following 
we will implicitly assume that the configurations y G y under consideration are smooth except 
in the segment y~x, ■ ■ ■ ,Vk for some K < N. We divide Q into an atomistic subdomain fi at 
such that yj G f2 at for all j G {— K, . . . , K} and the continuum domain Q ch = Q\Q at . In O cb we 
will use the Cauchy-Born approximation on a cell- by-cell basis. In f2 at we will use the atomistic 
model with Dirichlet boundary conditions as discussed in Section [3| 

This basic setting gives rise to a variety of possibilities including the precise choice of c?Sl at 
and the boundary conditions imposed on the atomistic subproblem. Both will in general depend 
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y-K-i y-K vk vk+i 



Figure 5.1. An illustration of the first a/c method. In Q at = (ol(^), clr(u)) 
the atomistic problem is solved with the Dirichlet boundary conditions g*{y)- 
Outside il at the Cauchy-Born approximation is used in all cells Qj. 



on the configuration y. Our main objective for £ ac is the existence of a weak formulation in 
the sense that 



D£ ac (y)-u = J a* c (x)Vu(x)dx, 



where u € S#(y) is a piecewise linear interpolant of u G U and cr ac is a stress function to be 
determined. If this weak formulation can be obtained, the consistency analysis reduces to error 



estimates on fields, as already seen in Lemma 4.4 



Throughout this section, G iTi(J2) denotes the solution of the original minimization prob- 
lem (2.1) for a given configuration y Gy. 

5.1. An a/c method with optimal boundary conditions. We place the boundary a of the 
atomistic subproblem halfway between the interface atoms, that is a = a(y) = o_r(?/)] T , 
where 

/ \ y-K-i + y-K ( s vk + vk+i 
o-lkv) = 2 ' aR ( y > = 2 ' 

Let O at = aii(y)) and il cb = S7\r2 at . We write the a/c energy £ &c {y) as the sum of a 

continuum and an atomistic part 

£ ac (y) = £?(y) + £*\y), (5.1) 

which are introduced below. 

Due to the choice of a(y) there are two half cells, (y-K-i, a^(y)) and (an(y), y^+i), in the 
continuum region il cb (see Figure 5.1). Since the cell problems are symmetric, the Cauchy- 
Born energies of these half cells are given by an d h^K+l(y)j respectively. Hence, the 

continuum contribution to the energy £ ac is defined by 

-K-l N 

£* h (y) = £ £t{y) + ¥ ch K (y) + \£ c K + i{y)+ £ £?(*)• (5-2) 

j=-JV+l i=A'+2 

The coordinates of the atoms in the atomistic region O at are represented by 

2/at = {V-K, ■■■ ,yK) T - 



For the definition of £ at (y) we consider the minimization problem (3.2 ) on the atomistic domain 
il at subject to the Dirichlet boundary conditions g*(y) = [g*i{y) <£r(2/)] T - In correspondence 
with Remark [4] and Section 3.3 they are given by 

*/. x 1 lL{y)+Tj R (y) 1 rj L (y) + -y R (y) 

9L(y) = -, — 7— » 9R(y) = -, — 7— > 

1 — T 1 + T 1 — T 1 + T 
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where r = e~ e a( - y) , and 7l,7_r are denned in (3.15). The energy contribution from the 
atomistic subproblem is thus given by 

£?{V) = £a(y),g*(y){y a t) = ~ ^ Uaiy)^, Vat) ■ V G tf 1 ^), plan- = 5%)}, 



where I a (y) is defined as in (3.3). We denote the solution of this optimization problem by 
4>Xt E -ff 1 (O at ). It satisfies the boundary-value problem 

-e 2 A0; t + mVat = Py m ft at , 

4>lt\d^ = 9*{y)- 

From a computational point of view g* (y) is also a convenient choice since this is equivalent to 



homogeneous Neumann boundary conditions. In Section 3.3 we deduced a clear interpretation 
of the effect of this choice of boundary data: besides the interaction among themselves, the 
atoms in f2 at interact with mirror atoms outside fi at . This is closely related to the geometric 
reconstruction idea for classical potentials [T7], 14"]. 



In analogy to (2.3) we search for minimizers of the total potential energy 

Ef(y)=£™(y) + {f,y) £ (5.3) 

in y, where / E U~ 1,2 represents an external force. Formally, a minimizer y ac satisfies the 
following Euler-Lagrange equation in U~ 1,2 : 

DEf(y) = D8 &c (y) + / = 0. 

Throughout the remainder of this article we assume that the atomistic domain f2 at is large 
compared with e, that is Aa 3> e and hence terms of order 0(t) are exponentially small. To 
keep the notation more compact we will not give precise estimates of r-dependent terms arising 
from the atomistic domain explicitly but include an 0(r) where necessary. 

5.2. Consistency. In order to study the consistency properties of the a/c energy £ ac (y) from 



(5.1) we first need to calculate its derivative. Having established weak formulations for the 
derivatives of £ , £ ch , as well as £ a ,g, we will prove that the a/c energy £ ac admits a similar 
reformulation of D£ ac (y)-u. For this we have to take into account that both the boundary of the 
atomistic domain f2 at and the boundary conditions depend on y. The necessary preparations 
were carried out in Section [3j 

Lemma 4.1. Let y E y satisfy miny' > Furthermore, let u £ U be a test vector and 
u E S#(y) an interpolant of u; then, 



f ( a ch (x) if x E O cb 

D£ ac (y).u= / a™(x)Vu(x)dx, where <(*) = I) \ J 

Jn I a y,* \ x ) y x G ^ 



(5.4) 



and cfA^x) is given by (2.7) with (ft = </>* t . 



Proof. 1 . Continuum Contribution. From Section [4] we already have the equality 

D£f{y)-u= [ a y h j (x)Vu(x)dx, 

j E {-N, -K — 1} U {K + 2, . . . , N}. For the contribution ^^(y) from the half cell 
(y-K-i, cbL(y)) we make use of the symmetry of the cell problems. Since Vu\q_ k is constant, 
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o-l{d) is the midpoint of Q-k = (y-K-i,y-K), and o y _ K is symmetric in Q-k, we deduce 
that 

\D£%{y)-u = - f a£_ K (x)Vu(x)dx = P " a y h _ K (x)Vu(x) dx. 
1 J Q-k Jy-K-i 

We treat \£f? +l (y) analogously. Hence, 

D£?(y)-u = [ a y h (x)Vu(x)dx, 



cr^j(x) if X G Qj. 



where cr y h {x) 

2. Atomistic Contribution. To calculate the derivative D£^(y) we use the chain rule and 
the derivatives that were provided in Section[3} Applying Proposition 3.1 (with ht = (u-k-i + 
u_k)/2, hji = (% + Uk + i)/2 because of D y a(y) ■ u = a(u)), we get 

D£f(y)-u = D y £ a{y)ig , {y) (y at )-u at + D a £ a{y) ^ (y) (y at )-D y a(y)-u 

(5.5) 

a^(x)Vu{x)dx, 



where the stress is given by (2.7) with 4> = 0* t and w a t = (u-k, ■ ■ ■ , uk) G M 2ii " +1 is 
the section of u corresponding to the atoms in the atomistic region. Note that the choice of 
boundary conditions implies D g £ a ^ g ,^(y at ) = 0; cf. RemarkEJ □ 



Remark 5. The weak form (5.4) of the derivative D£ AC already implies that there are no 
ghost forces for homogeneous deformations y. If the atoms are equidistant, then g* L {y) = </>(ol) 
and g\{y) = (t>{cLR) and thus also (p* t = (p in f2 at . Moreover, it is clear that tjjV' = <p for all j. 



Hence, we obtain that cr y c (x) = cr y (x) for all x £ £1, which implies that D£ ac {y) = D£(y) = 

for all y = FX G y representing homogeneous deformations (i.e., that the method exhibits no 
ghost forces). □ 

Absence of ghost forces does not immediately imply consistency of the a/c method, but has 
to be shown separately. This we do next. Because of the structure of the weak formulation 



(5.4), the analysis boils down to estimating the errors between the field 4> coming from the 
original atomistic model and the fields respectively, 0* t . 

Theorem 5.1. Let y G y be such that miny' > so > so/ then, for all u with interpolants 
u G S#{y), 

\{D£{y)-D£™{y))-u\ < C{e\\y" W^+r) \\Vu\y, 
where C = C(sq) and the weighted l 2 w -norm is defined by 



with weights Wj := max (l, e" 



II "Il2 

\\y Wfi 

W,SQ 

-msodist(j,{-K,K})} 



I // 1 

£ z^j=-N w j\yj\ 



(5.6) 



Proof. Using the weak formulation (5.4) of D£ a,c (y) we obtain 



{D y £(y)-D y £ ac (y)).u\ 



o- 



y (x) - al c (x))Vu(x)dx 



< Wa. 



y-Vy IIl 2 (o) 



N 

I £ e\\rr 

j=-N 



y 



iVnl 



qci|2 



L 2 



y nL°°(Qj) 



1/2 



iVttl 



L 2 



(5.7) 
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For Qj belonging to the continuum region Lemma 4.4 and Lemma 4.2 imply 



P-Qj-nj+n-l]) 



ne 



-mnsQ 



n=l 

oo 



<CeJ2\\y"\\iH\J-n, j+ n-i])n 3/2 e- mnso 



n=i 

oo 



< rVl \^ll?i"ll 2 „ p -mns 
S Cf^ 2^" y We 2 ({j-n,j+n-l]f 



71=1 



1/2 



(5.8) 



where we have employed the Cauchy-Schwarz inequality twice and used the fact that the series 
YlnLi n 3 e~ mns ° is convergent. 

Summing over all cells belonging to the continuum region and interchanging the order of 
summation we obtain 

oo 

,qci|2 .. ,- , .:>, V" \ " ll„."l|2 

je{~N,...,N} n=l 
\{-K+l,...,K} 

N 



e K 



rr-'-n- < rV 3 IUf"ll 2 -mns 



je{-N,...,N} 

\{-K+l,...,K} 



k=-N 



where 



E E • 

je{-N,...,N} n=l,...,oo 
\{-K+l,...,K} ke[j-n,j+n-l] 



-mnsQ 



This is a geometric series from which we can factor out e ms odist(A;,{ K,K}) ^ an( j ij ence we obtain 
w' k < Cwk, which gives 

N 

£ £ K - <W 2 l-( Qj) < Ce 3 £ (5-9) 
j${-K+l,...,K} k=-N 

To compute the consistency error of the weak form in the atomistic region, we need to bound 
the difference \\a y — ay C 1 1 Loc ^ ^ = \\a y — || i(X >(Q .) f° r Qj C ^ at . Using the same arguments 
as in the proof of Lemma 4.4 we obtain 



I at II ^ /~i 

\ a y a V,*\\L°°(Qj) - U 



£tlU«(Qj) + e ll V </> - V^|| L ao(Q )). 



Lemma 4.1 implies 

|| at || 

\\ a y a y,*\\L° 



e mm xSQj d a (x) 



m < C(|0(a L ) - <£(y)| + |0(a fl ) - g R (y)\)t 
Next, we recall from Remark [4] that 9* R = 1* R + C( r )) which is given by (cf. Remark [4| 

1R(V) = [ p T f{x)G £ (a R -x)dx + 0(t), 

where p V y^{z) = Yljez ^e( z ~~ V^) an< ^ ^ refl ^ s a reflected and periodized extension of {yj)f = _ K - 
Hence, we obtain 



{a R )-g R {y)\ = \<p{a R )- lL {y)\+0{T) 
1 



< 



2me 



{ P y(z)-pf(z))e- r -f^dz 



+ 0(t). 
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Minor modifications of the proofs of Lemma 4.2 and Lemma |4,1| yield 

„ oo 

\<Ka R ) - g R (y)\ < - E \yf ~ y^'" {min{yT ^- aR) + 0(r) 

£ j=K+l 

oo 

< CeJ2\\y"\\m[K-n+i,K+n])ne- mnso + 0(t). (5.10) 

71=1 

An analogous result holds for \4>{ai^) — g* L {y)\- It is now straightforward to see that the consis- 
tency error committed in the atomistic region can be bounded above in the same way as the 
consistency error committed in the continuum region (in fact it is dominated by (5.9). This 
completes the proof. □ 

5.3. Stability. The special choice g*(y) of boundary conditions for the atomistic subproblem 
allows for an elementary stability analysis of <? ac (y) that draws from the ideas we used in 
Section 13.31 We recall that 
1 

4me 



+ e -T^ +z - 2a ^)p y (z) dzdx + 0(t). 

The next result addresses the differentiability of and j R . We show that the derivatives 
satisfy certain bounds. 

Lemma 5.2. Let y G £l 2K+1 satisfy yi + \ — y%> for all i G {— K + 1, . . . , K}, a R — yx > 

eso/2, and y-K — cll > £?q/2. Then, 72 (y) is twice continuously differentiable with respect to 
y and a and there exists C(mmin y') (independent of e) such that 



\D 7L (y,a)-(u,h)\ < C(mmmy')^ U - K £ Hl J + £ K) 2 )^ 



k=-K+l 



\D 2 lL (y,a).[(u,h),(u,h)}\<C(mmmy')(( U - K + £ ^ 

V V - / U— IS 



k=-K+l 

„2 



for all u and h G M . Analogous bounds hold for j R (y, a). 

Proof. The proof is based on the observation that 

K „ K 



1 * — -v f m 1 \ ;/ m / \ « — -\ m 

7i(j,b) = - V / e-T^i-^di^e"! 1 "-"'' V e"7 

j=-K Jlla j=-K 
The rest of the proof is a straightforward computation; see [10t Lemma 5.3] for the details. □ 



The r-dependent terms in £^(y) = £ a (y)_ ! g*(y) (Vat) fr° m (3.24) only contain 72, (y) and j R (y), 
whose derivatives are bounded by Lemma |5.2| The derivatives of these r-dependent terms are 
therefore still of order 0(t) and will be neglected in the proof of the following result. 

Lemma 5.3. Let y G y satisfy miny' > %. Then, 

2 

D 2 S ac (y)-[u,u}> (^e- mm3Xy ' -0{Tfj\\u'\\% e Mu G hi. 
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Proof. We treat continuum and atomistic contributions independently and start with the for- 
mer. Lemma 14,51 states that 

D 2 £?(y)-[u,u] > ^e-^^Kl 2 



for all j = —N, . . . , N. Hence, the definition (5.2) of S% directly implies that 

2 2 /-K-l 1 N N 

rfSi\y).[uM>*— v,r ^e{ £ \u' 3 \ 2 + \{W_ K \ 2 + W K+l \ 2 ) + £ \u'A. 

^j=-N j=^+2 ^ 

Let us now turn to the atomistic part £% (y). From Section 3.3 we know that for the given 



choice of boundary conditions and a(y) we can write the energy of the atomistic part as 

K 

+ e -f {yK + 1+yK - x ' z) )S E (z- yj )dzdx 



Fit i — ^ 771 1 I 771/ s 

= i£_ u-^m-vA + & --^{yi+Vi-y-K-y-K-i) (5 11) 

i,j=—K m 

+ e -^(vK+yK + i-Vi-y^ + g sel{ + 0( r ), 

where the constant £ Be if accounts for the self-energies of the atoms {— K, . . . , K}. Differentiating 
twice and keeping only contributions from nearest neighbour interactions leads directly to 

D^y).[u,u}>e-"*y'^e( l -\u>_ K \ 2 + £ + \Wk + A ~ 0{r). 

\ i=-K+l ' 

Adding the lower bounds for D 2 £% h (y) • [u, u] and D 2 £^ t {y)- [u, u] we arrive at 

D 2 £™(y)-[u,u] = {D 2 £?(y)+D 2 £?(y))-[u,u] > (e"™^'^ - 0(r)) \\u'\\^, 

for all u £ U, as desired. □ 

5.4. Error Estimates. Combining the consistency and stability results we obtain the following 
error estimates. We note that the upper bound on the error depends on the smoothness of y in 
the continuum region, but that the dependence on y in the atomistic region decays exponentially 
with distance to the a/c interface. In realistic higher-dimensional models such an estimate would 
make it possible to allow defects in the atomistic region without affecting the error estimate. 

Theorem 5.4. Suppose that y £ sagmmEf and y ac £ arguim-E^ satisfy 

rainy', miny' ac > sq > and max y , max y ac < Sq < +oo. (5-12) 
There exist constants c and C = C(sq,So) such that, if Aa > clog(So), then 



\y'-y'^<C[s\\y"\\ Pwso +T). (5.13) 



Proof. From Lemma 5.3 it is clear that there exists a constant c such that, for Aa > c, we have 



D 2 £ &c (y)-[u,u]>^e- mSo \\u'\\ 2 l VueU, Vyey,y'<S . 



In particular, this holds for all y 6 conv{y,y ac }. Let cq = ^x~e 



my p -mS 
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Vk-j 



y K -i Vk a R (y) y K +i r] K+2 {y) 



VK+j+i[y) 



Figure 6.1. Illustration of the problem in the interval Qj = (yjc-j, 2a R (y) 
Vk-j) used to compute g R (y). 



Let u = y — y ac ; then we can choose y £ convjy, y ac } such that 

co|K||| < D 2 £™(y)-[u,u] = {D£* c (y) - D£^(yJ)[u]. 
Employing the consistency estimate of Theorem |5.1| we obtain the stated result. 



□ 



Remark 6. With some additional work it is possible to avoid assuming the existence of y ac , 
but deduce it from an inverse function theorem type argument |144 [TO] . □ 

6. Boundary Conditions From Cell Problems 



The boundary conditions g*(y) we imposed on the atomistic subproblem in Section 5.1 gave 
rise to a method without ghost forces, and whose analysis was relatively straightforward. The 



reasons for this is the clean weak formulation (5.4) of D£ ac and the convenient stability proper 



ties established in Lemma 5.3 We now investigate how this situation changes if computationally 
cheaper boundary conditions are chosen. The following construction may also provide a starting 
point for generalisations to higher dimensions. 

For example, a canonical choice, which requires no additional computational effort, is 

g L {y)=^- K \a L ) and g R (y) = ^ K+1 \a R ), (6.1) 

where we still assume az,(y) = \{y-K-\ + V-k) and a R (y) = \{vk + Vk+i)- In this case, we 
have the following result, which suggests that the additional error committed can be controlled. 



Lemma 6.1. Let m.m.y > sq > % and let gL/R be given by (6.1); then, 

\g(y)-g*{y)\<C{e l l 2 \\y"\\ ils +r). 



3.2) 



Proof. Without loss of generality we focus on g R only. Upon first estimating 
\g R (y) - gMv)\ < |V> ( * +1) M - <t>(a R )\ + \</>(a R ) - g* R (y)\, 



and then employing (5.10) and Lemma 4. 2 1 we obtain 



\g R (y)-g R (y)\<CeJ2\\y"\ 



^{[K-n+^K+n]) 



nc 



-mnsQ 



+ 0(t). 



n=l 



Using the same argument as in (5.8) to (5.9), we obtain the upper bound (6.2). 



Motivated by Lemma 6.1 we define a second a/c energy £ a,c (y) by 

£ ac (y)=£?(y)+£ at (y), 



where (y) is the same as in the method discussed in Section 5.1 (see (5.2)) and 



£ at (y) = £a(y),g(y)(y a t) 

= -inf {/«,(„)(¥>, !/ at ) : <p € H 1 ^ 1 ), tp\ao* = g(v)j- 



□ 



(6.3) 
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We denote the minimizer for given y by at 6 i^ 1 (S7 at ) . 

Before we embark on the analysis of this new method, we establish a useful auxiliary result. 



Lemma 6.2. Let mvay' > sq > Let g R {y) be defined by (6.1); then, it can be equivalently 
written as 



(6.4) 



In particular, g R is twice Frechet differentiable with respect to y, and there exists a constant 
C = C(sq) such that, for all y £ y with miny' > so > 

\ D y9n(y) ■ u \ < c Wk+i\ and \ D lgn(y) ■ [u,u]\ < C\u' K+l \ 2 Mu € U. 

Analogous results hold for gi(y)- 



Proof. Recall from (4.8) that 



9R(V) = ^ K+l \a R ) = ± £ [ S £ ( Z - y%£> 

Llll /id 



e e 



\a.R—z\ 



dz, 



where yj K+1 ^ denotes the periodic extension defined in (4.3). We use the identities 
V^K+k^ = VK + key' K+l and a R = y K + ^v'k+i 



to obtain 



9r(v) 



2m S L Si 



fceZ 



[z-Vk- key' K 



+i 



e e 



\yK+\ey' K+1 -z\ ^ 



Shifting the integration by (yx + \ £ u'k+i)^ we obtain (6.4). 

The bound on the first and second derivatives follows as in Lemma 5.2 the key observation 
being that 5 £ = 0(e^ 1 ) is balanced against the e preceding y'^+i i n its argument. □ 



6.1. Consistency. A crucial difference between the a/c energy (6.3) and the energy from 
Section 5.1 is that now the derivative of the atomistic energy with respect to the boundary 
conditions does not vanish. 



Since the continuum contribution to D£ ac (y)-u is the same as in Section 5.1 we only need 
to analyze D£ &t (y). Using the chain rule we obtain 

D£ at (y)-u = Dy at £ a{y): g {y) (y at )-u at + D a £ a{y)igiy) (y at )-a(u) 
+ D a S a(y), g (y)(y a t)-( D y9(y)-u). 



The same reasoning as in Section 5.1 gives for the first two terms on the right-hand side 



D y a ,t £ a.(y),g(y)(,y a ,t)- U a.t + D a£ a (y),g(y)(y a ,t)-a(u) 



(6.5) 



where a y (x) is given by ( |2.7[ ) with 4> = 4>ax- 

Next, we turn our attention to the term D g £ a r y \ g r y \(y a:t ) ■ (D y g(y) ■ u). We recall from 
Lemma 3.4 that (for Aa 3> e) 

D 9 £ a(y), g {y){y^) = -me[g L (y) - g* L {y), g R (y) - g* R (y)} + O(er). 
Combining this result with Lemma |6.2| and Lemma |6.1[ we obtain 

\ D 9 £ a(v),g(v)(v*k) ' ( D y9(v) - u )\< Ce(\g L - g* L \ + \g R - g* R \ + rj (\u'_ K \ 2 + |?4 +1 | 2 ) 1/2 

<c(e\\y"\\ ils +rV K|U, (6.6) 
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where C = C(miny'), and we have estimated £ 1 ' 2 t < r. Equipped with these estimates, we 
obtain the following consistency result. 

Lemma 6.3. Let y G y with miny' > sq > then, there exists a constant C = C(so) such 
that 



D£(y)-u-D£ a \y)-u 

where we have used the same notation as in Theorem\5.1 



<C{e\\y"h lso +r)\\Vu\\ L , 



VueU, 



Proof. From (|6.5|) and (|6.6) we obtain that 

< 



D£(y)-u-D£ at (y)-u 



£ 



_cb \\2 



.at||2 



j=-N,...,N 
H{-K+l,...,K} 



I a V a j,y \\l°°(Q j ) + \\ a y a V \\L2(a L ,a R ) 



1/2 



+ C(e\\y" 



\e 2 



|Vu|| L 2. 



The first group in the upper bound was already estimated in the proof of Theorem 5.1 and the 
second group, \\a y — o~^\ \L 2 ( a L ,a R ) can be treated analogously to the term \\a y — cr^\\L 2 (a L ,a R ) 



in the proof of Theorem 5.1 



□ 



6.2. Stability. We wish to compute a convenient lower bound on D 2 £' d,c {y) ■ [u, u] for some 
given y £ y with <io < sq < y' < Sq. Since the continuum part of the energy is the same as 
in the first method, we only address the stability of the atomistic subproblem with the given 
choice of boundary data. We write the second derivative of the energy £ at in the form 

D 2 £^(y) ■ [u, u] = D 2 £?(y)-[u, u] + (D 2 £ at (y) - D 2 £?(y)) ■ [u, u] 
and use the coercivity of D 2 £f J {y): we know from Lemma |5.3| that 

,2 /- K 



D 2 £*\y)-[u,u] > e ~"*y 



mjjL 



1 

2 £ V2 



/ |2 
U -K\ 



I '|2 | 1 / |2 

+ 2^ \ u *\ + n\ U K+l\ 



O(t) 



i=-K+l 

for all u 6 U; hence we are left to analyze the difference D 2 £ at (y) — D 2 £^ t (y). We will not 
show that this difference is small, but will only be able to bound it below by a controllable 
quantity. This is reminiscent of similar observations made in 



Lemma 6.4. Let y £ y such that rainy' > so > ?o; then there exists a constant C = C(sq) 
such that 

(D 2 £*\y) - D 2 £f(y)).[u,u] > -C{e l / 2 \\y"U w +r). 



Proof. The difference between the energies £ at (y) and £^{y) only consists of effects from the 



boundary conditions. We have, by (|3.18|), 

£ at (y)-£f(y) = 



Ia(y)(L(y),g(y),y) + Ia(y){Za{y),g*{y),y) = "y \9{v) ~ 9* (l/) \* + 0(st). 



As in Section 5.3, one can verify that the 0(r) term remains of that same order in the first and 
second derivatives. This implies that 

{D 2 £ at (y) -D 2 £?(y))-[u,u] = me{g(y) - g*(y)f \{D 2 g{y) - D 2 g*(y)) >,«]] 

|2 



+ 2me\ (Dg(y) - Dg*{y)) -u\ A + 0(t\\u'\\%) 



> me{g(y) - g* {y)f [{D 2 g{y) - D 2 g*(y)) ■ [u, u}] + 0(r\\u 



/||2 > 



3.7) 
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We now employ Lemma 6.2 to bound D 2 g(y), Lemma 5.2 to bound D 2 g* (up to another 0(t) 
error) , and Lemma |6.1| to bound g — g*, which yields 

(D 2 f a Hy)-^% at (y))-K«]>-G(£ 1/2 lly1le, S0 +^)ll u, ll|' 

where C = C(miny'). □ 



From Lemma [6. 4| and Lemma 5.3 we immediately obtain the following corollary, which states 
that, if So is moderate, y "smooth" in a neighbourhood of the interfaces a^/ji and in the 
continuum region, and if the atomistic region is sufficiently large, then D 2 £ ac {y) is stable. 

Corollary 6.5. Let y £ y satisfy mint/' > sq > ?o an d max?/ < Sq; then there exists a 
constant C = C(sq) such that 

,2 



D 2 £^(y).[u,u}>( 



m M e -mS 



V 2 



\u\\% 



Vu£U. 



Remark 7. The scaling e 1 ' 2 is due to the fact that the additional error committed is 
concentrated in a region of length e. □ 



6.3. Error Estimates. Repeating the proof of Theorem 5.4 verbatim, but replacing the con- 
sistency and stability estimates from Section [5] with those derived in Lemma 6.3 and Corollary 



6.5, we obtain the following error estimates for the modified a/c method. 



Theorem 7.1. Recall the notation introduced in Theorem 5.1. Suppose that y G argmini^ 
and y ac S argminE'f , where £ ac is defined in (6.3), satisfy 

miny', miny' ac > sq > and maxy',maxy' ac < Sq < +oo. (6.8) 

There exist constants c andC = C(sq,So) such that, if T+e l / 2 \\y"\\fi < ce~ mSo (in particular, 
K must be sufficiently large), then 

<c(e\\y"\\ £ 2 n t +t). (6.9) 



y 



7. Conclusions and Outlook 

We have presented a rigorous error analysis of an atomistic-to-continuum coupling method 
for a field-based interaction potential in one space dimension. The starting point for the design 
of coupling methods was a weak formulation of the forces arising from the atomistic model. 
This provided a natural connection point to the corresponding continuum model. We believe 
that the present work in a comparably simple setting addresses several important questions 
relevant for a/c coupling in the presence of fields, most prominently the dependence of the a/c 
methods on choice of the boundary and the boundary data for the interaction fields. 

For the two a/c methods we discussed we chose y-dependent boundaries a(y) of the atomistic 
subdomain $7 at . In other words we fixed the position of the boundary in the Lagrangian domain. 
This leads to convenient weak formulations of D£ ac (y). An obvious alternative (particularly 
relevant for higher dimensions) is the choice of y-independent a. We have not investigated this 
further, however, see [10] for some preliminary remarks. 

We also remark that we heavily utilized the one-dimensional setting in several places in the 
analysis. A generalisation both of the numerical methods and their analysis is therefore non- 
trivial. In particular, we can see no straightforward generalisation of the reflection boundary 
conditions g*{y ). A possible way forward would be to give an alternative analysis of the second 
method described in Section [6] that does not utilize these reflection techniques. 
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